Effects of Elliptical Polarization on Strong-Field Short-Pulse Double Ionization 
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We predict new end-of-pulse behavior in high-field atomic double ionization. Calculations of 
atomic electron trajectories in short intense laser pulses confirm our analysis of elliptical polarization. 
We exhibit a four-band structure in ion momentum distributions under various ellipticities, and 
predict that sequential and non-sequential double ionization can be cleanly distinguished under 
elliptical polarization. 

PACS numbers: 32.80.Rm, 32.60.+i 



When atoms or molecules are exposed to an intense 
(10 14 to 10 16 W/cm 2 ) femtosecond laser pulse, electrons 
that are released can collide under the laser force with the 
core. Many novel applications under the heading of at- 
tosecond physics are currently being vigorously pursued 
pp. Among the effects of importance is double ioniza- 
tion. This was discovered experimentally to exhibit ion 
yields 4-6 orders of magnitude greater than predicted by 
the usual sequential tunneling theory. A characteristic 
"knee" signature showing this enhancement has been ob- 
served with all rare gas atoms [2j E2 H] and also some 
molecules [5, 6 . The basis for the effect, i.e., for the ex- 
ceptionally strong electron pair correlation, is the subject 
of active research, both experimentally and theoretically. 

A three-step mechanism [7J illustrated in Fig. [I] has 
generally been accepted as descriptive of such double ion- 
ization events: (i) one electron is freed with zero kinetic 
energy by tunneling when the nuclear coulomb poten- 
tial has been tipped down by the strong laser field; (ii) 
the freed electron accelerates away in the laser field and 
is then driven back when the field reverses phase; and 
(iii) in returning, the laser-controlled electron acquires 
enough energy to carry away a second electron in a near- 
core e-e collision process. 

Many questions remain unsettled regarding strong- 
field double ionization, and one that is completely open 
concerns polarization. Almost all near-optical-frequency 
double ionization experiments have been carried out with 
linearly polarized light. The three-step picture suggests 
dramatically lower ion yield under elliptical or circular 
polarization, because the return trajectory is much less 
likely to encounter the atomic core, and this was quickly 
affirmed experimentally [8 . Double ionization has never- 
theless been reported even with circular or near-circular 
polarization in studies of atomic magnesium [9] and of 
several molecules [5]. 

The lack of cylindrical symmetry has up to now put 
the case of elliptically polarized fields beyond the range 
of essentially all quantum calculational approaches that 
are relevant to high-field double ionization [TOl [H] . How- 
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FIG. 1: The left panel shows the nature of the three-step rec- 
ollision process for highly correlated double ionization. The right 
panel shows the flat-top model for the field given in Eqn. (TIJ, here 
with <j> = and e = 0.5. The red and blue curves are the x and y 
components of the field. We take w = 0.0584 a.u., i.e., A = 780 nm. 
The 10-cycle pulse used in our calculations would have a duration 
of about 26 fs. We choose Eq = 0.131 a.u., corresponding to a peak 
(not average) intensity of / = 0.6 PW/cm 2 . 

ever, with the recent development of two- and three- 
dimensional approaches to high-field ionization problems 
via classical ensembles [12 , a viable theoretical avenue 
is open [13l HH [15] [TBI [17] . We report here surprising 
results from the first systematic theoretical exploration 
of double ionization with elliptically polarized light: 

E(t) = E f(t) [e x sin(wt + (f>) + e y e cos(wi + </>)], (1) 

as shown in Fig. [I] 

Since multiphoton experiments in this intensity regime 
examine recoil ion momentum distributions |10l 1181 1191 
12011111113 as well as ion yield, we have calculated numeri- 
cally the double ionization momentum distributions for a 
full range of ellipticities. A view of our classical ensemble 
results is provided in Fig. [2]below. They significantly ex- 
tend findings for non-zero ellipticity up to 0.5 reported by 
Shvetsov-Shilovski, et al. [23]. We also present an analyt- 
ical theory giving explanations for our numerical results 
and also for (a) an earlier speculative calculation [21] and 
(b) previously unexplained features seen in very high field 
double ionization experiments on neon |25j . We predict 
that the structures discovered should be observable with 
intensities and pulse durations readily available. 

Our classical ensemble method has been described 
previously [12] . and it has already contributed to in- 
terpretation of double ionization phenomena [26 . The 
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FIG. 2: End-of-pulse momentum distributions for doubly ionized ions with ellipticities from to 1, as indicated by the top label of each 
panel. 



Z = +2 charged ion core is fixed at the origin and 
in the energy of a two-electron bound system, E to t — 



\v\ + \v\ + V(fi,r 2 ), we take 



V(n,f 2 ) = -ZV sc (f u a) - ZV sc (r 2 ,a) + V sc (n - r 2 ,6), 

with V sc standing for the soft-core Coulomb potential of 
the model [271: 



V sc (r,c) = 



1 



V r 2 



(2) 



A many-pilot-atom method [28 is used to generate mi- 
crocanonical ground state ensembles. In the absence of 
the laser field, the energy E to t of each 2e member of the 
ensemble is set to be -1.3 a.u. We expect that this value 
is close enough to the ground state energy of both Kr 
and Xe to be useful in guiding experiments [29 . The 
parameter a can be regarded as the way the model ac- 
counts for the atomic core electrons. We set a — 1.77, 
to prevent autoionization, and b — 0.1, to allow strong 
electron pair correlation. Prior work with the model [13 
has shown that out-of-plane effects are negligible in first 
approximation at the intensities used experimentally, so 
we restrict the electrons to motion in the x-y polarization 
plane for simplicity. The end-of-pulse momentum distri- 
butions of doubly ionized ions, which are calculated as 
the sum of the momenta of the two ionized electrons, are 
shown in Fig. [2] 

Let us look at each graph in detail. For linear polariza- 
tion there is one continuous region in the momentum do- 
main, which lies along the p^-axis. As ellipticity increases 
to 0.2, this region grows into two parts: the right-upper 
part that lies mainly in the first quadrant and the left- 
bottom part that lies mainly in the third quadrant. As 
ellipticity increases to 0.3, each of the two regions fur- 
ther grows into two new parts and a four-band structure 
emerges. At ellipticity 0.5, the four-band structure can 



be clearly seen. Besides, a new elliptical structure be- 
comes obvious, part of which appears as two "bridges" 
connecting the four-band structure. This elliptical struc- 
ture grows and in the meantime the four-band structure 
decreases as ellipticity increases. At e = 1 the distribu- 
tion is circularly symmetric. 

To determine the origin of these momentum distri- 
butions, the history of each two-electron trajectory was 
traced back. Trajectory back analysis [3T] provides great 
insight into physical processes and is only possible with a 
classical approach. Depending on whether there are close 
recollisions or not, events are classified as non-sequential 
double ionization (NSDI) or sequential double ionization 
(SDI), as shown in Fig. U This fi gure is the first to show 
distinct patterns of ionization in which SDI and NSDI arc 
separated cleanly by ellipticity. Our observation that in 




FIG. 3: Momentum distribution of DI (top) can be further di- 
vided into momentum distribution of NSDI (left) and momentum 
distribution of SDI (right). Ellipticity used here is 0.5. 
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FIG. 4: Momentum distribution of SDI in x-direction (left) and 
in y-direction (right). The P x distribution shows a single peak 
centered at zero while the P y distribution shows four-peaks. Four 
vertical dashed lines are positions of peaks predicted by our analysis 
(see Eqn. (py), which clearly fits well with the numerical results. 



linear polarization 10% of the single ionizations are con- 
verted to SDI events appears consistent with 1-D reports 
from the phase space perspective [17] . 

One can see in Fig. [3] that SDI is responsible for the 
four-band structure and NSDI for the elliptical structure. 
It is helpful to project the 2D momentum distribution of 
SDI onto the x- and y-directions. The SDI distribution 
along the x-direction is a single peak structure centered 
at zero while the distribution along the y-direction has 
four peaks centered at approximately ±0.5 a.u. and ±1.5 
a.u., as shown in Fig. [3] This four-peak character has 
probably already been recorded without notice at least 
once in very high field sequential ionization of neon under 
nearly circular polarization (see the study by Maharjan, 
et al. [25]). 

Before beginning an analysis of the predicted four- 
peak structure one important feature needs to be pointed 
out. For single or sequential double ionization, elec- 
trons are almost always ionized [3D] along the x-direction 
(plus or minus) when the field is near maximum because 
the peak field in the x-direction is higher than that in 
the y-direction, due to ellipticity. At those x-maximum 
times, the field in the y-direction is zero. As a con- 
sequence, the electron momentum distribution shows a 
double peak structure centered at ± 0.5 a.u. along the 
x-direction while a single peak is centered at zero along 
the y-direction, as indicated by Fig. [5] 

After the electrons are ionized, the reasonable approx- 
imation can be made that the Coulomb attraction be- 
tween the electron and the ion core can be neglected and 
the only force that acts on the freed electron is the laser 
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FIG. 5: Momentum distributions of electrons (first or second) 
immediately after ionization along the x- and y- directions. It 
is a double peak structure centered at ± 0.5 a.u. along the x- 
direction while a single peak structure is centered at zero along the 
y-direction. 




x x 

FIG. 6: In-phase and out-of-phase (NZ and Z) trajectories are 
shown for clarity. 

force. Then the change of momenta of the two electrons 
due to drift motion with the laser field can be derived 
analytically. We let t\ and ti be the times when the first 
and second electrons are ionized. Then the total change 
of momentum of the two electrons due to drift motion 
along the y-direction is 

AP y = Apiy + Ap 2y 

= [wti smtojti + cf>) + 47T sin(wt2 ± <M — cos tj>] 

E s 

= Voti smhoti ± (f>) ± 47r sin((jjii + </>) — cos </>] . 

As we explain below, sin(wti + cj>) ~ 1 and the final cos <f> 
should be dropped, in which case we find four preferred 
values, to be compared with Fig. [4j 

AP y w 0.089 x (2-7T ± 4tt) am(wti + <f>). (3) 

Now let us explain the origin of the terms in our analy- 
sis in detail, and evaluate them numerically. Because we 
are in a regime where a significant amount of SDI can be 
achieved, it is obvious that the first electron will be ion- 
ized before the plateau stage of the field, i.e., during the 
turning-on stage, since the first electron is much easier 
to ionize than the second one. This point is confirmed 
by our trajectory calculations. It turns out that most of 
the first electrons are ionized around the beginning of the 
second cycle, so iot\ is in the near neighborhood of 2tt and 
cos cj> averages to zero. The difference between ti and t\ 
could be any number of half cycles - 2ir/u>, 3tt/uj, Att/uj 
and so on, as indicated by arrows in Fig. [1] In the odd 
half-cycle cases, when t^ — ti — (2n±l)7r/o;, the two elec- 
trons leave the ion core along opposite directions (along 
+x and —x). The result has been called [32] the "Z" sce- 
nario, for approximately zero momentum transfer, or out 
of phase double ionization. In the even half-cycle cases, 
ti — t\ = 2mt/uj, the two electrons leave the ion core 
along the same direction (either +x or —x) and we label 
this "NZ" for substantially non-zero momentum transfer, 
or in phase double ionization. 

The interesting y-components of the end-of-pulse mo- 
menta of the two electrons are explained as follows. For 
in phase SDI, the y-components of the end-of-pulse mo- 
menta of the two electrons most probably also have the 
same sign, as shown in Fig. [6] Hence the magnitude of 
the y-component of the end-of-pulse momentum of the 
ion is large, which corresponds to the outer two peaks 
of the four-peak structure in Fig. [4] For out of phase 
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SDI the y-components of the end-of-pulse momenta of 
the two electrons most probably also have opposite sign, 
hence the magnitude of the y-component of the end-of- 
pulsc momentum of the ion is small, which corresponds 
to the inner two peaks of the four-peak structure. On 
the assumption that the electrons are ionized exactly 
when the field in the x-direction is maximum, we have 
sin(cjii + 0) = ±1 and sm(ujt2 + 4>) — sin(wti + 2n7T+0) = 
sin(wii+0) (in phase) or sin(wi2 + <f>) = sin(u>ii + (2n + 
1)tt + 4>) = — sin(oj£i + <j))(out of phase). The peak 
positions predicted in Eq.(|3j are ±0.56, ±1.68, close to 
matching the ensemble average values shown in Fig. [4j 

In summary, we have predicted a new end-of-pulse 
momentum distribution for sequential double ionization, 
and have shown that elliptical polarization provides a 
new control mechanism for recollision physics in high- 
field ionization. Our trajectory calculations are made 
in the spirit of the standard recollision picture and are 



analytically feasible because of our special trapezoidal 
pulse shape. The numerical distributions are only pos- 
sible to obtain classically at this time. This is because 
no workable non-perturbative quantum approach to el- 
liptical high-field short-pulse effects is known. The good 
agreement achieved provides additional support for the 
interpretation of high-field ionization processes as exam- 
ples of electron physics first, and atomic physics secon- 
darily, a view that encourages classical interpretations 
and calculations of even more complex high-field effects. 
Finally we remark that the same four bands with the 
same separation are obtained numerically from short 
smooth Gaussian pulses as well. 
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